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Abstract. We report the discovery of a new quadruply imaged quasar surrounded by an optical Einstein ring candidate. Spectra 
of the different components of IRXS J 1 13155.4-123155 reveal a source at z = 0.658. Up to now, this object is the closest known 
gravitationally lensed quasar. The lensing galaxy is clearly detected. Its redshift is measured to be z = 0.295. Additionally, the 
total V magnitude of the system has varied by 0.3 mag between two epochs separated by 33 weeks. The measured relative 
astrometry of the lensed images is best fitted with an SIS model plus shear. This modeling suggests very high magnification of 
the source (up to 50 for the total magnification) and predicts flux ratios between the lensed images significantly different from 
what is actually observed. This suggests that the lensed images may be affected by a combination of micro or milli-lensing and 
dust extinction effects. 
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1. Introduction 

The peculiar and complex morphology of the source 
IRXS Jl 13155.4-123155 (hereafter J1131GL) has been 
serendipitously unveiled during polarimetric imaging of a sam- 
ple of radio quasars carried out in May 2002 at ESO, La 
Silla. These observations are reported in Section |2l together 
with additional optical imaging obtained in December 2002. 
Astrometry and photometry of the gravitational lens system are 
also described. In Section|3 we present spectroscopic observa- 
tions of the source, the lens and a nearby companion (here- 
after Jl 131b). Sectionl^is devoted to a simple lens model and 
Section|5lsummarizes why this new gravitational lens is a par- 
ticularly interesting one. We have adopted throughout the paper 
Ho = 65 km s"' Mpc"', = 0.3 and Aq = 0.7. 

2. Direct imaging, astrometry and photometry 

Direct imaging of J1131GL has been obtained at two differ- 
ent epochs with respectively EFOSC-2 at the 3.6 m telescope 
and EMMI-Red at the 3.5 m New Technology Telescope (NTT) 
at the La Silla observatory. On May 2, 2002, we resolved the 
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four components of J1131GL on a set of 2x4 polarimetric 
images (corresponding to 4 different orientations of the Half 
Wave Plate and 2x150 s integration time per orientation) taken 
through a combined V-band + WoUaston prism. The average 
seeing measured on the frames is 1.1" and the pixel size is 
0.158". Additional V and R images of J1131GL have also been 
obtained under poor seeing conditions (1.6") on December 18, 
2002. The coadded exposure time amounts to 480 s in V and 
960 s in R. The pixel size is 0.166". 



2.1. Image analysis 

By means of the GENERAL image decomposition program 
developed by Remy et al. (|1997 ), we have fitted the EFOSC 
images of J1131GL (see Fig.^) using 4 free adjustable PSF 
components, altogether with a galactic light profile for the lens. 
Due to significant distortions of the PSF across the field, we 
have finally obtained the best results using a purely analytical 
2D Moffat PSF 

Fitting a (non seeing-convolved) de Vaucouleurs or expo- 
nential disk profile allows to roughly remove the low frequency 
signal of the lensing galaxy and to unveil a ring-like struc- 
ture at a level higher than 3cr above the noise. This feature 
is spectacularly confirmed by the deconvolution of the images 
using the MCS code (Magain et al.'1998'(: a non uniform ring 
passing through the 4 point-like components is seen on the de- 
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Fig. 1. (a) V direct image of Jl 13 IGL obtained with EFOSC-2 (4x150 s integration time). The 4 lensed images and the deflecting 
galaxy are identified on this CCD frame, (b) Deconvolved image using the MCS method (see text) and (c) Idem but with the 
quasar images removed. A conspicuous Einstein ring candidate is seen on the deconvolved images. 

Table 1. Left: Relative positions between the different lensed components (B, C, D) and the lensing galaxy (G) with respect to 
component A as deducted from the best seeing EFOSC-2 frame using the GENERAL (GEN) and MCS codes (see text). Right: 
Relative magnitude of the B, C, D lensed images with respect to A deducted with the MCS method at two different epochs. The 
last two columns give respectively the amplification //exp (sign - parity) predicted by the SlS+y model and the corresponding 
relative magnitude with respect to A (Section|3. 
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convolved images (Figs.H'iC). Additionally, we measured that 
the integrated signal from the gravitational lens system (A-D 
and G) was not significantly polarized (i.e. total polarization ~ 
0.2+0.08%). GENERAL and MCS have also been applied to 
the EMMl frames but despite of a more stable PSF across the 
field, the poorer seeing prevented us to reach better results than 
with EFOSC-2. 



rors on the photometry and on the astrometry are probably not 
negligible for D and G due to their small separation and their 
equally low S/N. 



2.3. Absolute photometry and variability 



2.2. Relative astrometry and photometry 

Astrometry of the lensed components (B-D) and of the lens- 
ing galaxy (G) relatively to A has been derived using the 
GENERAL and MCS codes. With the MCS algorithm, the po- 
sition of the lens G has been measured by deconvolving a point- 
like source at the lens position. The results listed in Table [2 
were derived for the best seeing EFOSC frame. They are in 
statistical agreement with the positions retrieved for the sec- 
ond epoch. The image deconvolved by the MCS code is a sum 
of analytical point sources and of a diffuse numerical decon- 
volved background. Due to this, the ring and the lensing galaxy 
are better taken into account with the MCS code rather than 
with GENERAL. Consequently, the resulting photometry of 
the point-like components obtained with MCS is more reliable. 
The results obtained for the two epochs are reported in Tabled 
The 1 cr errors listed in Tableware formal errors on the fit for the 
results derived with GENERAL. In the case of the MCS code, 
they reflect the dispersion of the results when changing the ini- 
tial conditions and deconvolution parameters (i.e. smoothing, 
position, flux and background step in the;^^^ fit). Systematic er- 



The two polarimetric standards HD 155197 and HD 298383 
have been observed with the same setting (i.e. V band -i- 
Wollaston) during the photometric night on May 2, 2002. These 
stars allowed us to calculate a zeropoint of 25.85+0.03 mag. 
Consequently, the integrated V magnitude of the system is esti- 
mated to be 16.63 and Va = 17.97+0.09 (deducted with MCS). 
The brightest parts of the ring have a surface brightness V ~ 
23.9 mag arcsec"^. The estimated mean surface brightness of 
the galaxy inside a 4" radius is V ~ 22.7 mag arcsec"^ and its 
integrated magnitude inside the same radius is V=18.4. Since 
the conditions were not photometric during the December ob- 
servations, we performed differential photometry of the inte- 
grated system J1131GL with respect to various objects in the 
field. Due to the polarimetric nature of the EFOSC frames (i.e. 
the field is splitted in non contiguous bands), only 4 objects 
could be used for this purpose; and it was necessary to also use 
the 10 s acquisition frame. We found that the integrated flux 
of the system (A-D and G) was brighter by 0.29+0.04 mag in 
December 2002. Since the relative photometry between the 4 
components is quite similar at both epochs (see Table in- 
trinsic variability is very likely responsible for this diff'erence. 
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Fig. 2. (a) Integrated spectrum of the three bright components 
A, B, C. (b) Spectrum of A (2x900 s) and D (900 s). Emission 
hnes of a quasar at z = 0.658 are easily identified, (c) Spectrum 
of the lensing galaxy (900 s). The typical absorption lines of an 
elliptical galaxy at z = 0.295 are identified. The flux scales are 
arbitrary. The small horizontal lines above the wavelength axis 
on (c) indicate telluric absorption lines. 



3. Low resolution spectroscopy with ElVIIVII 

Two sets of low resolution spectra (2" slit) have been obtained 
using the new CCD on the EMMI-Red arm at the NTT; (1) at 
the end of May 2002, we obtained a 900 s unresolved spec- 
trum (PA ~ 15°, from North to East) of the three bright com- 
ponents A, B, C with the CCD in the 1x1 bin mode and (2) on 
January 25, 2003, we obtained 2x900 s spectra wifli flie CCD 
in the 2x2 bin mode (corresponding to 3.58 A pixel"') and the 
slit passing through A, D, G and J1131b, located at 25" East 
of J1131GL. Standard bias subtraction, flatfielding and spec- 
trum extraction procedures were used. Wavelength calibrations 
of the spectra were done with an He-Ar lamp. Due to the ab- 
sence of a spectrophotometric standard in January 2003, we 
used the spectrophotometric standards LTT 1788 and LTT 2415 
observed with a 5" slit in December 2002 to correct these data. 

The integrated spectrum of A, B, C obtained in May 2002 
(Fig.Et) shows a continuum and emission lines typical of a red- 
shifted quasar. We could separate these three components us- 
ing the deconvolution algorithm based on a Maximum Entropy 
Method developed by Khmil and Surdej (2002i. The sepa- 
rated spectra of A, B, C obtained by this algorithm were found 
to be very similar with slight possible changes in the con- 
tinuum slope but they also remained highly correlated due to 
the poor seeing of 1.7". The second set of spectra obtained in 
January 2003 under an average seeing of 0.66" have completed 
our knowledge of this system. Due to the spatial sampling of 
0.332"pixel the spectra of A, G and D slightly overlapped. 
Nevertheless a careful choice of the apertures enabled us to ex- 
tract the spectrum of A on both frames and the spectra of G and 
D on the second one. An average spectrum of A is shown in 
Fig-Eb- We identify on this spectrum typical emission lines of 
a quasar at a redshift z=0. 658+0. 001. This redshift was deter- 
mined by fitting a Gaussian on the Mgll A 2798 A and [OIII] AA 



4959, 5007 A emission lines. Unfortunately, the 2" slit slightly 
scatters the light from B and C (located less than 1.2" away 
from A) and consequently the spectrum of A is contaminated 
on each exposure by these two components. 

Because of the overlapping and of the similar brightness 
of G and D their spectra are mutually contaminated. We thus 
subtracted one spectrum from the other after adequate scaling 
in order to visually minimize the contamination (i.e. to remove 
quasar emission lines from G and the 4000 A break from A). 
This handicraft process suggests a relative reciprocal contam- 
ination smaller than 10%. The resulting spectra are shown in 
Figs. EJi-c. Component D shows similar emission lines as A, 
but has a flatter slope, not generated by decontamination. The 
spectrum of G is typical of an elliptical galaxy (Kennicutt ll992> 
showing absorption lines redshifted at z = 0.295+0.002 (e.g. 
Call K&H AA 3933, 3968 A, G band A 4304 A, etc). 

The spectrum of Jl 131b shows Balmer absorption lines and 
a continuum typical of an A type star contaminated by a back- 
ground galaxy lying at 2" from this object. This star was previ- 
ously identified by Bauer et al. J2000> as the most likely source 
of the X-ray (ROSAT) and radio (NVSS) emission present in 
this field. They also reported a redshift z=0.654 for this object. 
We firmly reject their identification, mistakenly matched with 
the quasar we have observed here. 

4. A simple lens model 

The relative positions of the quasar images and of the lens- 
ing galaxy with respect to image A (see Table Q have been 
fitted using two simple lens models: the Singular Isothermal 
ElHpsoid (SIE, Kormann et al. 1994. Kassiola & Kovner[T993} 
and the Singular Isothermal Sphere plus an external shear 
(SlS+y). The best fit is obtained with the SIS+y model and 
yields the following results: the angular Einstein radius Be - 
1.819+0.006", the shear y = 0.123+0.003 and the shear posi- 
tion angle cp - 14.84+0.11° (from North to East). The uncer- 
tainties on the parameters come from the fit of the model on 
1000 Monte-Carlo synthetic observations compatible with the 
observed errors. The direction orthogonal to the shear axis (p 
does not point towards any bright object close to the lens. The 
reduced (for 3 degrees of freedom) is significantly smaller 
for the SIS+r (Xy - 19) than for the SIE lens model (xli^ - 
203) because the latter cannot reproduce correctly the observed 
lens position. If one does not fit the lens galaxy position, the re- 
duced is comparable for the SIE and the SIS+y but the SIE 
model predicts a lens position 0.4" away from the observed 
one. Finally, using the formulae of Witt et al. (|2000) for sin- 
gular isothermal lens models with shear, we predict the time 
delays between B and the other multiple images: Tbc - 0.01 d, 
Tba - 0.87 d and Tbd - 96 d. The leading image is B and the 
time delay sequence is thus BCAD. 

As can be seen from Fig.|3j this system is a long axis quad 
with a source located very close to a cusp (at about 0.1"). This 
allows very high magnification of the host galaxy at a few hun- 
dreds of pc from the center. The latter is the likely origin of 
the optical ring-like structure joining the lensed point-like im- 
ages. The two minima images of the arrival time surface, B 
and C (with positive parity), should have about half the flux of 
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Fig. 3. Results of the SIS+y model expressed in a system 
of coordinates centered on the lensing galaxy G. Left panel: 
the source plane position of the source (x), relative to G, is 
(-0.411", -0.1084")- This corresponds to a closest distance to 
a fold caustic of 0.016" (115 pc in the source plane). Right 
panel: the contours (calculated from the model) corresponding 
respectively to 110, 740, 3700 and 5200 pc from the central 
engine in the source plane are overplotted on the deconvolved 
EFOSC frame. 

the central one (the saddle point A), as demonstrated by Mao 
nWlS and by Schneider and Weiss (TWl. The observed B/A 
magnification ratio (see Table [!) does not follow this generic 
prediction: B is even the brightest image. Since the C/A ratio 
is "normal", and since the time delays are very short between 
A, B, C, we suspect B to be affected by micro/milli-lensing (D 
could also be). Even if the V-band and R-band flux ratios are 
quite similar (see Tabled, we cannot exclude that dust extinc- 
tion also plays a role in the observed flux ratios. The unlensed 
absolute magnitude of the source is Mb = -22.7 if we conser- 
vatively take - 10 and B - V - 0.2. Thus, the source is, 
strictly speaking, an AGN/Seyfert 1 . 

5. Discussion and conclusions 

Direct imaging and long slit spectroscopy of 
IRXS Jl 13 155.4-123 155 presented in this letter have 
enabled us to show that this object is a quadruply imaged 
quasar {V,„i e [17.5,20.1]; AO ~ 1.2") at redshift z = 
0.658±0.001 lensed by an efliptical galaxy at z = 0.295±0.002. 
The MCS method has made possible to enhance the signature 
of an Einstein ring candidate barely seen on the direct images. 
The simple modeling by an SIS +y shows that the source is 
located very close to a cusp allowing the host galaxy to cross 
the caustic and to generate the ring seen on the deconvolved 
images. As it is observed in many quads, there is a discrepancy 
between the observed image flux ratios (especially between 
the saddle point A and the minimum B images) and the ones 
predicted by modeling. This may reflect the necessity to use 
a lens model involving a small percentage of substructures 
(Schechter and Wambsganss 120021 references therein). More 
data on J1131GL are necessary before drawing any definite 
interpretation of the observed discrepancy. Nevertheless, the 
explanation of the flux ratios in this system should probably 
involve micro/milli-lensing and/or dust extinction. Only 
individual spectra obtained simultaneously from UV to NIR 



(and taken at time intervals equal to the time-delays) for 
each component will enable one to disentangle between these 
effects. Finally, we have shown that the integrated flux has 
varied by 0.3 mag between May and December 2002. 

The source J1131b located at 25" East of J1131GL was 
wrongly identified by Bauer et al. ( I2000t as the most likely 
optical counterpart of the X-ray ( Voges et al. [T999j and radio 
(Condon et al. 1998 ) emission present in this region. Our iden- 
tification of this object as a hot star suggests that Jl 131GL (al- 
ready cataloged with the RXS notation) is the true source of the 
X-ray and radio emission. Note that the discovery of a grav- 
itationally lensed system in a multi-wavelength survey is not 
surprising due to the large expected multi-band magnification 
bias (Borgeest et al. 1991 Wvithe et al. l2003t . 

This first set of data suggests the necessity of good spa- 
tial resolution and high signal to noise ratio multi-wavelength 
imaging (from radio to X-ray) in order to use the rare charac- 
teristics of this system as many observational constraints for an 
accurate modeling. In this framework, the shape of the Einstein 
ring is an invaluable asset to determine independently the shape 
of the lens potential and of the unlensed source (Kochanek et 
al. l2001i) . Moreover, a precise lens inversion of the lens equa- 
tion (e.g. Warren and Dve 120031 should enable one to draw a 
unique multi-wavelength picture of the source and to retrieve 
information on the source at angular scales inaccessible with 
present day and even future instrumentation. 

This bright system brings together rare properties (i.e. 
quad, bright optical Einstein ring, small redshift, high ampli- 
fication), nearly unique among the known gravitational lens 
systems. These features make IRXS Jl 13155.4-123155 a very 
promising astrophysical laboratory for future investigations, in- 
cluding the possibility for an independent determination of the 
Hubble parameter Ho based on time delay measurements. 
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